HZPP-9807 
Dec. 10, 1998 



SCALING EXPONENTS AND FLUCTUATION STRENGTH 
IN HIGH ENERGY COLLISIONS 



Liu Lianshou, Fu Jinghua and Wu Yuanfang 

Institute of Particle Physics, Huazhong Normal University, Wuhan 430079 China 



ABSTRACT 

The information on dynamical fluctuations that can be extracted 
from the anomalous scaling observed recently in hadron-hadron col- 
lision experiments is discussed in some detail. A parameter "effec- 
tive fluctuation strength" is proposed to estimate the strength of dy- 
namical fluctuations. The method for extracting its value from the 
experimentally observed quantities is given. Some examples for the 
application of this method to real experimental data are presented. 
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The JACEE experiment in 1983 W and afterwards some accelerator experiments' 2,3 ' 
indicated that there are dynamical fluctuations beyond the usual statistical ones in 
high energy multiparticle final states. These experimental phenomena had caused high 
interest in studying the anomalous scaling in high energy collisions' 4 ' . 

Recently, the expected anomalous scaling has been observed successfully in both the 
data of 7r + p and K + p collisions at 250 GeV/c from NA22' 5 ' and the data of pp collisions 
at 400 GeV/c' 6 ' from NA27, under the assumption that the dynamical fluctuations are 
anisotropic in longitudinal-transverse planes' 7 '. 

It is now natural to ask the question: What is the underlying physics of these 
experimental findings, and/or what are the properties of the dynamical fluctuations 
that lead to the experimentally observed anomalous scaling. This question is the 
subject of a long term study. In this letter we want to try a limited task in this 
direction, i.e. to study the relation between the observed anomalous scaling exponents 
and the strength of dynamical fluctuations within the framwork of a simple model - 
the random cascading a-model. 

Firstly, let us recall that what is observed directly in experiments is the anomalous 
scaling of normalized factorial moments F q , which under the assumption of Poisson 
or Bernoulli type of statistical fluctuations are equal to the normalized probability 
moments' 8 ' : 

1 M (n q ) 

Here, a phase space region A is divided into M sub-cellsQ and p m is the probability for 
finding particles in the mth sub-cell. The anomalous scaling of F q : 

F,(M)ocM*« (M -> oo) (2) 

with non- vanishing indices <p q , called intermittency (IM) indices, is an evidence for the 
existence of dynamical fluctuations' 8 '. 

On the other hand, in studying the fractal property of the system, one usually uses 
the un-normalized moments' 9 ' 

M 

c q =Y,(p m q )> (3) 

m=l 

which has the anomalous scaling property 

C q {M)(xM- {q - 1)D " (M -> oo). (4) 

The indices D q , called multifractal Renyi dimensions' 9 ', are related to the IM indices 

4> q as 

D 1 =1 ~ -ftl- (5) 

1 Throughout this paper we use the convention popular in the study of intermittency and fractality 
in high energy physics, i.e. let M in Eq's. (1) - (4) be the number of sub-cells in the division of 
the phase space region A of topological dimension Dt- This means, in particular, that for a normal 
geometrical object the fractal dimension defined via Eq.(4) is always equal to unity. On the contrary, 
in the usual definition of fractal dimension' 9 ' : C q (l) oc Z^ -1 )- ? (/ — » 0), I is the 'edge' (1-dimensional 
scale) of the hypercube (sub-cell), and concequently for a normal geometrical object D q = Dt- 



The question in front of us is as the following: When we obtain successfully a strict 
scaling of normalized factorial moments through a proper way of phase space division, 
as for example in the cases of Ref.'s [5] and [6], we are ready to get the slopes of 
the \nF q ~lnM plots, which gives the IM indices and Renyi dimensions. Then, how 
to extract a characteristic quantity describing the strength of dynamical fluctuations 
from the values of these indices and dimensions. This is the problem we are about to 
discuss. 

It should be noticed that an anomalous power law and a fractal property of a 
system both result from the dynamical fluctuations in this system, and therefore the 
IM indices <p q and multifractal Renyi dimensions D q are related to the strength of 
dynamical fluctuations. They have, however, their own physical meaning and cannot 
be taken directly as a measure of the dynamical-fluctuation strength. For example' 9 ', 
the first-order Renyi dimension Di, or equivalently the information dimension Dj, 
determines the scaling property of the number of boxes containing the dominant part 
of information; the second-order Renyi dimension D2, sometimes called correlation 
dimension /i, measures the scaling properties of two particle correlations, etc. These 
scaling properties come from a common origin — the dynamical fluctuations. All of 
them are related to these fluctuations but none of them can serve as an appropriate 
quantity for describing the strength of these fluctuations directly. 

In order to show the relation between the strength of dynamical fluctuations and 
the values of IM indices or Renyi dimensions, let us consider a simple example — the 
random cascading a-model' 8 ' 10 '. This model describes each multiparticle event as a 
series of steps, in which the initial phase space region A is repeatedly divided into 
A = 2 parts. After v steps we get M = 2 V sub-cells of size S = A/M. At each step v 
the normalized particle density is obtained in each of the two parts by multiplication 
of the normalized density in the step v — 1 by a particular value of the random variable 
u) U j v1 where j u is the position of a sub-cell at the uth step (1< j u < 2 U ). 

The elementary fluctuation probability 00 can be chosen in various ways' 8 ' 10 !. The 
simplest way that provides a characteristic parameter for describing the fluctuation 
strength is to choose it as' 10 ! 



in which r is a random number distributed uniformly in the interval [—1,1], j is an 
integer (1 < j < 2 Jy_1 ), a is a positive constant taking value in the range [0,1], 



w v ,2j-i = -(l+ar) 



\{ l -<*r), 



(6) 



< a < 1 . 



(7) 



The value of a determines the possible region of u, 



1 — a 
2 



< u < 



l + a 



2 



Let us note that u defines the way how particles are distributed from step to step 
between the two pieces of a given cell, i.e. it characterizes the strength of multiplicity 



fluctuations in cell-division, and a determines the width of the possible values of tv, 
therefore, a is the characteristic quantity describing the strength of dynamical fluctua- 
tions in this version of random cascading model. In the following, we will analyse the 
relation between the strength parameter a and the multifractal dimensions D q in this 
model. 

In the random cascading a-model 

v 

In the limit of large z/, the distribution of random variable ( = ln^(«) approaches a 

i=i 

Gaussian' 8 ' 

p(()d( = (2ttz/)- 1 /V- 1 exp[-(C - vV) 2 /2va 2 ]d( 

with 

a 2 = J P{uj){\a{uj)-V) 2 duj , V = j P(uj)\n(uj)duj. (9) 
Explicit calculation of multifractal dimensions in this limit gives 

D , = 1 - ^\ (10) 

A characteristic feature of the Gaussian approximation is the proportionality of 
1 — D q and q as can be seen from Eq.(10): 

1 ~ D q _ 1 j2 
q 2\n2 a ' 

Here and in the following we take A = 2 for simplicity. In Fig.l is shown the relation 
between 1 — D q and q for a = 0.1 to 0.5. It can be seen from the figure that when a is 
small there is a fairly good linear relation between 1 — D q and q, so that the Gaussian 
approximation is sufficiently good in these cases. 

To get a relation between D q and a, we calculate the variance a 2 , appearing in 
eq.(10), of the random variable lnw 

a 2 = (\n 2 uo)-(\nuo) 2 = \a 2 + \a A + --- . (11) 
Under linear approximation 

a 2 » l -a 2 . (12) 

How a is related to a is shown in Fig. 2. In this figure, the full circles represent the 
result without linear approximation, the dashed line indicates the linear-approximation. 
When a is not very large (a < 0.5 say), the two results are nearly equal. Therefore, 
both Gaussian and linear approximations can be used when a < 0.5. This region of a 
is sufficient for the limited range available in actual experiments, cf. Table I. 



Substituting Eq.(12) into Eq.(lO), we get 



1 2 



D q = 1 — — — -qa 



61n2 

or 

a = 



61n2, . , . 

-£)„). (13) 

It can be seen from Eq.(13) that, in the random cascading a-model, with small a, the 
strength parameter a of dynamical fluctuations is related to the multifractal dimensions 
D q by a very simple relation. Using Eq.(13) we can get an approximate value of the 
fluctuation parameter a as long as the multifractal dimension of any order q is known. 
In high energy experiments, the second-order multifractal dimension D2 is m ost easy to 
obtain. For this reason, the r.h.s. of Eq.(13) for q = 2, ^3 ln2(l — D 2 ) ~ ^2(1 — D 2 ), 
can be taken as a characte ristic quan tity for the strength of dynamical fluctuations. 

In Fig. 3 the values of \J2(1 — D 2 ), denoted by a c s (cf. Eq.(14)), are ploted against 
the model parameter a varying fr om to 1. The dashed line corresponds to a e s = a. 
From the figure we can see that \J2(1 — D 2 ) has almost equal value with a, especially 

when a is not very large. This means that \J2(1 — D 2 ) represents the value of a 
fairly well. Thus, we have been successful in obtaining an estimation of the fluctuation 
strength a in terms of the second-order fractal dimension D 2 in the framework of the 
random cascading a-model. 

The above results are obtained from a special model. In the general case, when 
the underlying dynamics is unclear, it is hard to define the "strength of dynamical 
fluctuations" strictly. In these cases, we can make use of the results obtained from the 
random cascading a-model as an estimation for this strength. 

Thus, for an arbitrary process that has anomalous scaling property, we define an 
effective fluctuation strength 



a eS = y/2{l-D 2 ) = yfefa. (14) 

as an estimation of the strength of the dynamical fluctuations taking place in this 
process. Its physical meaning is: 



The effective fluctuation strength a e g of an arbitrary process is the fluctuation 
strength of a random cascading a-model with elementary partition number A = 2 
that can give the same value of second-order IM index (fi 2 (within the Gaussian 
and linear approximation) as this process. 

In getting the last equality of Eq.(14) the relation (5) between the IM indices and 
Renyi dimensions has been used. 

Using the effective fluctuation strength defined above, we are now able to compare 
the strength of dynamical fluctuations in different collision processes. In the following 
we will give some examples. Before doing that, a question has still to be considered. 



In real experiments the dynamical fluctuations exist in higher-dimension t 11 ' and are 
usually anisotropic^ 5- ^ with a particular value of Hurst exponent Hn± . How can we use 
the definition (14) of effective fluctuation strength, which depends on a one-dimensional 
a-model with elementary partition number A = 2, to these cases? 

In answering this question, let us note that in case of anisotropic dynamical fluc- 
tuation (self-affine fractal), the partition numbers in different phase space directions 
cannot be simultaneously equal to integer values. The method of factorial moment 
analysis with non-integer partition^ 12 ' 13 ! has to be used and the resulting F| D for ar- 
bitrary value of M 3D all lie on a same straight line. Therefore, we can freely choose 
M 3D = 2 U ,v = 1,2, . . .. The anomalous scaling property of such a 3-D fractal is equiv- 
alent to that of a 1-D fractal with elementary partition number A = 2. Therefore, our 
definition (14) for effective fluctuation strength is applicable also to this case. 

Now, let us turn to the examples for the application of effective fluctuation strength 
to real experimental data. For this purpose we have to choose those data that possess 
good anmalous scaling property. 

For hadron-hadron collisions, the presently available 3-D data that have good scal- 
ing property are the self-affinely analysed data for 7r + p and K + p collisions at 250 GeV/c 
from NA22' 5 ' with Hurst exponents H ypt = H yip = 0.475 , H ptV = 1. The second-order 
IM index is obtained as [13] : = 0.061 ± 0.010. 

Another example of hadron-hadron collisions is the pp collisions at 400 GeV/c 
from NA27 [6] . In this case, only 2-D (r], (p) data are available due to lack of momentum 
measurement. The Hurst exponent in the (rj, (p) plane is found to be H vip = 0.74. 
Fitting the results to a straight line gives 2 . = 0.051 ± 0.004. 

As an example of e + e~ collisions we take the data from DELPHI^ 14 '. After omitting 

the first point to eliminate the influence of momentum conservation^ 15 ' , a good fit to 

a straight line comes out0, cf. Fig. 4. The second-order IM index is then obtained as 
03D = omQ ± o 005 

The resulting effective fluctuation strengths a e s for these three cases are listed in 
the last column of Table I. Its physical meaning is that, the anomalous scaling of the 
second-order factorial moments in these 3 experiments can be produced by random 
cascading processes having fluctuation strength a ~ 0.349, 0.319, 0.446, respectively. 



Table I The Hurst exponents, IM indices, second order Renyi dimensions and 
effective fluctuation strengths for the data from 3 experiments 



Experiment 


Hurst exponent 


02 


D 2 


«efT 


NA22 (3D) 


H y± = 0.475, H ptv = 1 


0.061±0.010 


0.939±0.010 


0.349±0.028 


NA27 (2D) 


H w = 0.74 


0.051±0.004 


0.949±0.004 


0.319±0.014 


DELPHI (3D) 


H= 1 


0.099±0.005 


0.901±0.005 


0.446±0.012 



Let us notice that all the values of 02 and D 2 listed in Table I are near by the 
boundaries of their allowed range, i.e. for 2 and 1 for D 2 . This would give us 

2 The reason why the e + e~ data have good scaling property already for the Hurst exponent H = 1 
will be discussed elsewhere. 



an impression that the strength of dynamical fluctuations in these experiments are all 
marginal. This is, however, wrong. From the last column of Table I we can see that the 
values of a e g for hadron-hadron collisions are approximately equal to 1/3 while that for 
e + e" collisions is close to 1/2. Since the allowed range of the parameter a characterizing 
dynamical fluctuations is [0, 1], cf. Eq.(7), the above results show that the dynamical 
fluctuations are nearly equal to one third of the maximum possible strength in hadron- 
hadron collisions and about half of the maximum possible strength in e + e~ case. This 
gives us, at least qualitatively, a feeling about the strength of dynamical fluctuations 
in these collision processes. 

Acknowledgements 

We are grateful to Zhang Yang and Liu Fuming for helpful discussions. This work 
is supported in part by the National Natural Science Fund of China. 



References 

[I] T. H. Burnett et al. (JACEE), Phys. Rev. Lett. 50 (1983) 2062. 

[2] P. Carlson (UA5), 4th Topical Workshop on pp collider physics, Bern, March 1983; 

G. J. Alner et al. (UA5), Phys. Rep. 154 (1987) 247. 
[3] M. Adamus et al. (NA22), Phys. Lett. B185 (1987) 200. 
[4] E.A. De Wolf, I.M. Dremin and W. Kittel, Phys. Rep. 270 (1996) 1. 
[5] N.M. Agababyan et al. (NA22), Phys. Lett. B382 (1996) 305. 
[6] S. Wang, Z. Wang and C. Wu, Phys. Lett. B410 (1997) 323. 
[7] Wu Yuanfang and Liu Lianshou, Phys. Rev. Lett 70 (1993) 3197. 
[8] A. Bialas and R. Peschanski, Nucl. Phys. B273 (1986) 703, B308 (1988) 857. 
[9] G. Paladin and A. Vulpiani, Phys. Rep. 156 (1987) 147. 

[10] Wu Yuanfang, Zhang Kunshi and Liu Lianshou, Chinese Science Bulletin 36 
(1991) 1077. 

[II] W. Ochs, Phys. Lett. B247 (1990) 101. 

[12] Liu Lianshou, Zhang Yang and Wu Yuanfang, Z. Phys. C69 (1996) 323. 

[13] N.M. Agababyan et al. (NA22), Phys. Lett. B431 (1998) 451. 

[14] P. Abreu (DELPHI), Nucl. Phys. B386 (1992) 471. 

[15] Liu Lianshou, Zhang Yang and Deng Yue, Z. Phys. C73 (1997) 535. 




b 



Sqrt(1/3) * a 



0.8 p- 

0.7 ~- 

0.6 ~- 

0.5 ~- 

0.4 ^ 

0.3 ~- 

0.2 r- 

0.1 - 



Q gj I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.10.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 2 Relation between the standard deviation a of random variable Incu 
and model parameter a. 
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Fig. 3 Relation between the effective fluctuation strength a e g and model parameter a. 
The dashed line corresponds to a cS = a 
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Fig. 4 The anomalous scaling of second order 3-D factorial moment of e + e collisions 
at 91 GeV (data taken from Ref.[14]). 




